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ABSTRACT 


The  Millstone  Hill  Thomson  (incoherent)  scatter  radar  system  has  been  operated 
routinely  since  1963  to  perform  a  synoptic  study  of  F-region  electron  densities, 
and  electron  and  ion  temperatures.  This  report  describes  system  changes  made 
in  1968  which  considerably  increased  the  accuracy  of  the  measurements  and  al¬ 
lowed  their  extension  to  higher  altitudes.  These  changes  have  also  made  it  pos¬ 
sible  to  measure  the  vertical  velocity  of  the  plasma  over  the  altitude  range  450 
to  900  km  to  an  accuracy  on  the  order  of  5  to  10  m/sec,  depending  upon  altitude 
and  time  of  day. 

Of  even  greater  significance,  complete  machine  reduction  of  the  results  is  now 
possible  so  that  considerable  savings  in  time  and  effort  have  been  secured  in  an¬ 
alyzing  the  data.  The  new  system  permits  all  the  radar  data  to  be  gathered  in 
the  digital  computer  in  real  time,  thereby  eliminating  the  need  for  post  real-time 
processing  of  magnetic -tape  recordings  of  the  signals.  Furthermore,  it  is  now 
possible  to  transmit  the  data  to  other  workers  in  computer-usable  form. 

This  report  describes  the  main  functions  of  the  computer  program  required  to 
analyze  the  measurements,  and  lists  the  times  of  all  measurements  made  with 
the  new  system  prior  to  1  January  1970.  Examples  of  these  results  are  presented 
and  discussed. 


Accepted  for  the  Air  Force 

Franklin  C.  Hudson 

Chief,  Lincoln  Laboratory  Office 
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INCOHERENT  SCATTER  MEASUREMENTS 
OF  F-REGION  DENSITY,  TEMPERATURES, 
AND  VERTICAL  VELOCITY  AT  MILLSTONE  HILL 


I.  INTRODUCTION 

A  vertically  directed  incoherent  scatter  radar  system  operating  at  68-cm  wavelength 
(440  MHz)  was  brought  into  operation  at  Millstone  Hill,  Westford,  Massachusetts  (42.6°N,  71. 5°  W) 
in  1963.  This  radar  system  has  since  been  in  routine  use  securing  measurements  of  F-region 
electron  density,  and  electron  and  ion  temperatures.  The  results  obtained  during  1963  through 
1965  have  been  published  in  earlier  technical  reports  as  well  as  in  a  number  of  journal  articles 
listed  in  Tabic  I. 


TABLE  1 

PUBLICATIONS  CONCERNING  THE  MILLSTONE  HILL  UHF 
(68-CM  WAVELENGTH)  THOMSON  SCATTER  RESULTS 

Year 

Months  Covered 

Publicotian 

1963 

March,  July,  August,  September 
December,  April,  July,  November 

February  1963  to  January  1964 

Ref.  2 

Ref.  3 

Technical  Repart  374, 

Lincoln  Laboratory,  M.l.T. 

(22  January  1965),  DDC  AD-616607 

1964 

April,  July,  November 

January  through  December 

Ref.  4 

Technical  Report  430, 

Lincoln  Laboratory,  M.l.T. 

(15  November  1967),  DDC  AD-668436 

1965 

Jonuory,  April,  August 

January  through  December 

Ref.  5 

Technical  Repart  474, 

Lincoln  Laboratory,  M.l.T. 

(8  December  1969) 

The  determination  of  the  electron  and  ion  temperatures  at  a  given  altitude  depends  upon  a 
spectrum  analysis  of  the  signals  returned  from  that  altitude.  Prior  to  1968,  we  accomplished 
this  with  a  spectrum  analyzer  that  could  examine  only  a  single  height  at  a  time.  Accordingly, 
throughout  the  period  January  1965  to  July  1968,  we  used  a  magnetic-tape  recorder  to  store  the 
received  signals  for  later  analysis.  In  this  way,  a  considerable  increase  in  the  time  resolution 
could  be  achieved,  but  it  was  necessary  to  spend  a  great  deal  of  time  in  processing  the  magnetic 
tapes.  This  scheme  had  other  serious  drawbacks  —  most  notably,  the  separation  of  the  results 
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for  density  and  temperature  made  it  necessary  to  perform  the  final  stages  of  data  reduction  by 
hand;  this  proved  to  be  most  time  consuming  and  unattractive  work. 

Beginning  in  1966,  we  seriously  considered  the  construction  of  a  spectrum  analyzer  that 
could  examine  the  signals  from  all  altitudes  of  interest  at  the  actual  time  of  receipt,  i.e.,  in  real 
time.  We  also  considered  it  desirable  that  the  integration  of  these  signal  samples  be  drift-free, 
and  yield  results  that  could  be  stored  in  the  digital  computer  along  with  the  measurements  from 
which  the  electron  density  is  determined.  In  1967,  we  began  the  construction  of  a  system  capable 
of  performing  these  tasks,  and  were  able  to  place  it  into  operation  in  July  1968. 

In  Sec.  II,  the  new  spectrum  analyzer  is  described,  as  well  as  the  other  changes  that  have 
been  made  to  the  equipment  since  it  was  first  brought  into  operation.  Section  III  describes  the 
operating  procedure  for  the  radar  equipment  developed  since  July  1968  that  provides  semiauto¬ 
matic  control  of  the  system;  also  given  in  this  section  are  tables  of  times  for  which  measure¬ 
ments  have  been  made  prior  to  1  January  1970.  The  functions  of  the  computer  program  that  has 
evolved  to  handle  the  data  analysis  are  described  in  Sec.  IV.  Development  of  this  program  re¬ 
quired  considerable  effort,  and  took  at  least  as  long  as  the  development  of  the  new  spectrum 
analyzer.  Section  V  presents  and  discusses  briefly  some  results  produced  by  the  new  system. 

II.  EQUIPMENT 

A.  System  Parameters 

The  main  features  of  the  incoherent  scatter  radar  system  have  been  described  previously.^'^ 
Although  numerous  changes  and  improvements  have  been  made  to  the  equipment  since  it  was 
brought  into  operation,  the  system  parameters  have  not  changed  appreciably  —  with  the  exception 
of  the  noise  temperature  of  the  radar  receiver,  which  has  been  lowered.  Table  II  lists  the  param¬ 
eters  of  the  system  as  of  1  January  1970. 


TABLE  II 

RADAR  SYSTEM  PARAMETERS  AS  OF  1  JANUARY  1970 

Antenna  diameter 

220  feet 

Antenna  collecting  area 

-1600  m2 

Polarization 

Transmitted 

Circular 

Received 

Opposite  sense  circular 

Peak  transmitter  power 

3. 0  MW  maximum 

Average  transmitter  power* 

120  kW  maximum 

Pulse  length  (See  Table  IV)* 

2.0  msec  maximum 

System  temperature* 

~180°K 

Receiver  pre-detector  bandwidtht 

40  kHz 

Receiver  post-detector  time  constantt 

50  psec  far  100-psec  pulses 

200  psec  for  0.5-,  1.0-, 

and  2.0-msec  pulses 

*  Monitored  continuously. 

t  Far  the  power  profile  channel  (Sec.  Il-E). 
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As  noted  in  Table  II,  many  of  the  system  parameters  are  monitored  continuously,  and  the 
antenna  collecting  area  can  be  determined  from  observations  of  the  intensity  of  the  radio  source 
Cygnus  A  which  transits  the  beam  daily.  However,  the  measurements  of  electron  density,  elec¬ 
tron  and  ion  temperatures,  or  drift  velocity  do  not  depend  upon  absolute  measurements  of  echo 
intensity,  and  the  monitoring  of  the  various  parameters  of  the  radar  is  undertaken  solely  to  en¬ 
sure  that  it  is  functioning  properly. 

B.  Antenna 

The  radar  antenna  consists  of  a  fixed  vertically  pointing  220-foot -diameter  parabola  with  a 
focal  length  of  100  feet  (see  Ref.  1).  Since  its  completion,  the  antenna  has  been  surveyed  annually 
and  adjustments  made  to  maintain  a  true  parabolic  contour.  Some  of  the  concrete  post  supports 
have  settled,  requiring  that  they  be  replaced. 

Apart  from  this  and  other  routine  maintenance,  the  only  changes  made  to  the  antenna  system 
were  undertaken  in  an  attempt  to  minimize  ground-clutter  echoes.  In  1967,  we  placed  a  large 
conical  skirt  around  the  feed  horn  in  order  to  cut  down  the  intensity  of  any  signals  transmitted 
toward  the  horizon.  With  this  done,  we  considered  scattering  from  the  tripod  legs  of  the  feed 
support  system  as  the  largest  remaining  source  of  unwanted  radiation  at  low  angles.  Accordingly, 
we  placed  absorbing  netting  on  the  insides  of  these  legs  to  eliminate  direct  reflections.  Forward 
scattering  around  the  legs  was  still  possible,  however,  and  to  reduce  this  it  was  necessary  to 
cut  down  the  energy  incident  on  them.  We  accomplished  this  by  constructing  an  extension  for  the 
feed  horn  that  increases  the  taper  of  the  illumination  over  the  dish  from  10  to  20  dB.  This  ex¬ 
tension  also  serves  to  cut  down  the  radiation  at  large  angles  by  virtue  of  a  quarter-wave  choke 
that  surrounds  it.  No  loss  of  antenna  gain  was  observed  when  the  extension  was  installed  —  prob¬ 
ably  because  the  horn  efficiency  was  raised  by  the  presence  of  the  choke. 

Although  these  changes  have  not  eliminated  ground-clutter  echoes,  they  have  lowered  their 
intensity,  and  it  is  no  longer  necessary  to  blank  the  receiver  out  to  ranges  (~150km)  where  the 
last  echoes  occur.  Indeed,  the  clutter  echoes  do  not  now  saturate  the  receiver,  so  that  it  is  pos¬ 
sible  to  attempt  to  remove  them  in  the  data  processing,  we  are  developing  programs  that  attempt 
to  do  this. 

C.  Receiver 

The  receiver  system  is  a  straightforward  multiple  superheterodyne,  having  intermediate 
frequencies  of  30  and  2.0  MHz,  that  is  preceded  by  a  parametric  amplifier.  All  the  local- 
oscillator  frequencies  are  derived  from  the  station  master  oscillator  to  achieve  good  frequency 
stability.  This  is  required  for  measurement  of  the  vertical  drift  velocity  of  the  plasma  (Sec.  1V-D) . 

The  original  receiver  preamplifier  was  an  electron-beam  parametric  amplifier  (Adler  tube) 
that  was  synchronously  pumped.  This  device  gave  extremely  reliable  service,  and  was  essen¬ 
tially  immune  to  large  amounts  of  RF  leakage  from  the  transmitter.  Synchronous  pumping  (at 
880  MHz)  yielded  the  lowest  system  noise,  but  folded  the  radar  signals  about  their  center  fre¬ 
quency  and  made  it  impossible  to  search  for  frequency  shifts  produced  by  drifts.  In  1967,  the 
Adler  tube  was  replaced  by  a  regenerative  diode-parametric  amplifier  operating  at  ambient  tem¬ 
perature  and  employing  a  klystron  oscillator  to  supply  pump  power.  This  amplifier  yields  a 
lower  system  temperature  (Table  11)  as  well  as  eliminating  the  fold-over  problem.  The  vacuum 
tube  receiver  has  been  replaced  by  a  completely  "solid-state"  receiver  in  which  gating  is  accom¬ 
plished  with  commercially  available  diode  gates.  With  the  possible  exception  of  the  parametric 
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amplifier,  the  new  system  has  greater  gain  stability  than  the  old  system.  On  the  other  hand,  the 
parametric  amplifier  has  been  found  to  be  somewhat  temperature  sensitive.  Day-to-night  tem¬ 
perature  changes  in  the  hut  (under  the  antenna)  in  which  the  amplifier  is  located  required  that  it 
be  retuned  several  times  during  a  24-hour  period.  By  improving  the  thermal  insulation  and  air 
conditioning  the  hut,  this  problem  has  largely  been  remedied. 

A  predetector  bandwidth  of  40  kHz  is  employed  in  the  receiver,  and  this  is  set  by  a  filter 
placed  in  the  2 -MHz  IF  channel.  Unlike  the  earlier  system,  this  same  bandwidth  is  employed 
for  all  pulse  lengths,  and  only  the  post -detector  time  constant  is  varied  (Table  II). 

The  receiver  system  temperature  is  continuously  monitored  by  including  in  the  sweep  inte¬ 
gration  process  a  calibration  pulse  consisting  of  a  known  amount  of  RF  noise  power  inserted  at 
the  receiver  terminals  via  a  directional  coupler.  This  noise  pulse  also  serves  to  establish  the 
absolute  power  in  the  echoes. 

D.  Transmitter 

The  transmitter  has  not  been  modified  significantly  since  it  was  built.  However,  we  recently 
replaced  the  vacuum  tube  exciter  with  one  employing  solid-state  components.  This  exciter  is  a 
reverse  superheterodyne  employing  the  same  intermediate  frequencies  as  the  receiver.  Mutual 
interference  is  kept  to  a  minimum  by  extensive  screening,  and  by  maintaining  all  signals  at  low 
power  levels  until  the  final  frequency  of  440  MHz  is  reached. 

Although  the  transmitter  is  rated  at  150-kW  average  power  output,  we  have  found  that  it  is 
not  usually  possible  to  exceed  130-kW  output  before  the  maximum  permissible  beam  current  is 
reached.  At  the  low  duty  cycle  employed  in  many  of  the  measurements,  it  is  possible  to  achieve 
a  peak  power  of  3  MW,  i.e.,  exceeding  the  2.5-MW  rating  of  the  transmitter.  However,  no  ad¬ 
verse  effects  have  been  encountered  operating  at  this  peak  power. 

E.  Power  Profile  Channel 

Measurements  of  the  distribution  of  echo  power  with  range  are  made  by  means  of  sweep  in- 
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tegration  using  a  digital  computer.  This  measurement  is  termed  the  power  profile  and  leads  to 
a  determination  of  the  electron  density  distribution  with  altitude. 

Power  profile  measurements  are  carried  out  by  sampling  the  receiver  output  repeatedly  at 
equal  intervals  along  the  timebase,  and  summing  together  all  samples  taken  at  the  same  position 
on  successive  sweeps.  For  this  purpose,  the  2-MIIz  IF  signals  are  converted  in  frequency  to 
200  kHz  and  then  applied  to  a  linear  rectifier  and  sampled. 

The  rectifier  is  designed  to  be  linear  over  the  range  0.1  to  10  volts,  and  is  followed  by  an 
DC  low-pass  filter.  Two  low-pass  filters  are  available  having  time  constants  of  50  and  200  psec. 
Selection  between  them  is  carried  out  by  a  relay.  The  filtered  rectified  signal  is  applied  to  an 
analog-to-digital  converter,  possessing  1024  levels,  that  operates  over  a  range  of  10  volts.  For 
observations  made  with  0.1 -msec  pulses,  samples  are  taken  at  50-psec  intervals,  and  the  50-psec 
low-pass  filter  is  used.  For  all  other  pulse  lengths,  the  sample  spacing  and  low -pass  filter  time 
constant  are  200  psec  (Sec.  I1I-A).  Samples  accepted  by  the  computer  are  immediately  squared 
before  being  summed.  In  this  way,  a  true  square-law  rectification  process  is  achieved  with  a 
dynamic  range  of  the  order  of  50  dB. 
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F.  New  Spectrum  Analyzer 


1.  Limitations  of  Original  System 

The  most  significant  change  made  to  the  equipment  has  been  the  construction  of  a  new  spec¬ 
trum  analyzer.  This  has  made  a  considerable  improvement  in  the  accuracy  of  the  results,  as 
well  as  rendering  all  the  data  analysis  open  to  machine  processing  as  noted  in  the  Introduction 
{Sec.  I).  The  original  spectrum  analyzer  and  later  modifications  have  been  described  else- 
where.  '  Basically,  these  systems  employed  a  bank  of  narrow -band  crystal  filters  to  analyze 
a  sample  of  the  received  signals  gated  from  the  radar  timebase.  The  duration  of  the  sample  was 
made  equal  to  the  length  T  of  the  transmitted  pulse,  and  the  delay  t  was  chosen  to  correspond 
to  a  height  of  interest.  The  filter  outputs  were  rectified  and  stored  in  analog  (Miller)  integrators. 

This  scheme  suffered  a  number  of  drawbacks.  It  was  possible  to  examine  only  a  single  de¬ 
lay  at  a  time,  and  thus  to  analyze  echoes  from  all  heights  of  interest  required  about  90  minutes 
observing  time.  In  order  to  reduce  this,  we  began  in  1965  to  record  the  IF  signals  on  magnetic 
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tape  for  later  playback  and  analysis.  This  proved  to  be  a  time-consuming  operation  since  ap¬ 
proximately  2  to  2\  weeks  were  needed  to  process  the  tapes  obtained  in  a  single  24-hour  period. 
Although  we  followed  the  best  possible  construction  practice  available  at  the  time,  the  analog  in¬ 
tegrators  were  subject  to  DC  drifts  that  served  to  introduce  "noise'1  into  the  measurements.  As 
a  result,  the  maximum  useful  length  of  an  integration  period  was  limited  to  about  5  minutes,  so 
that  the  minimum  signal -to-noise  ratio  for  which  useful  measurements  could  be  obtained  was 
about  0.1. 

A  further  limitation  was  imposed  by  the  filters  themselves,  which  did  not  have  equal  response 
above  and  below  the  peak  transmission.  As  a  result,  the  spectra  obtained  were  asymmetric  about 
the  radar  frequency,  and  it  was  necessary  to  average  the  two  sides  to  remove  the  effect."^  Finally, 
the  spectrum  analyzer  constituted  a  large  piece  of  electronic  equipment  requiring  constant  atten¬ 
tion  and  maintenance. 

2.  First  Plan  for  New  System 

In  searching  for  a  replacement  for  the  system  described  above,  we  sought  an  approach  that 
would  have  the  following  characteristics: 

(a)  Examine  all  heights  of  interest  simultaneously, 

(b)  Provide  accurate  drift -free  integration, 

(c)  Operate  in  real  time  (i.e.,  eliminate  analog  data  storage), 

(d)  Place  the  spectrum  results  in  the  computer  along  with  the  corresponding 
power  profile  results. 

A  scheme  which  achieves  these  goals  has  been  developed  at  the  Arecibo  Observatory  by  Perkins 
and  Wand£  In  this  scheme,  the  spectrum  analysis  of  the  signals  is  carried  out  by  forming  the 
autocorrelation  function  of  successive  samples  of  the  signals,  and  then  obtaining  the  Fourier 
transform.  Thus,  for  example,  to  obtain  the  spectra  of  signals  observed  with  a  T  =  0.5-msec 
pulse,  samples  of  the  received  signals  are  taken  at  20-psec  intervals  over  a  0.5-msec  interval 
of  the  timebase,  and  the  autocorrelation  function  of  these  samples  is  determined  in  a  digital  com¬ 
puter.  If  carried  out  in  full,  this  requires  the  computer  to  perform  3236  multiply  operations  for 
every  0.5-msec  interval  examined.  Savings  may  be  made,  however,  if  not  all  the  lagged  products 
are  formed.  In  the  limit,  only  26  multiply  operations  are  needed  to  form  one  sample  of  each  of 
the  possible  lagged  products  0,  20,  40,  .  .  .  ,  500  psec.  Perkins  and  Wand  call  this  the  "single** 
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sample"  method,  and  employ  this  approach  to  examine  altitudes  for  which  the  pre-detector 
signal-to-noise  ratio  is  expected  to  be  large,  i.e.,  where  the  samples  are  highly  correlated  so 
that  it  is  unnecessary  to  form  all  possible  products.  They  estimate^  that  the  loss  in  accuracy 
under  these  circumstances  is  equivalent  to  reducing  the  observing  time  by  about  a  factor  of  four. 
Thus,  by  reducing  the  number  of  multiply  operations  per  correlation  function  to  26,  it  is  possible 
to  examine  a  large  number  of  altitudes  simultaneously  at  the  loss  of  a  certain  amount  of  time 
resolution.  A  similar  tradeoff  may  be  had  by  adopting  the  "one-bit"  correlator  approach  as  dis¬ 
cussed  by  Farley7 

At  Millstone,  the  signal-to-noise  ratio  is  less  than  unPy  at  many  heights  of  interest.  Un¬ 
der  these  circumstances,  the  accuracy  of  the  correlation  method  depends  on  the  square  root  of 

Q 

the  number  of  lagged  products  that  are  formed,  so  that  the  single-sample  method  is  no  longer 
acceptable.  Accordingly,  we  decided  that  the  full  autocorrelation  function  would  be  computed. 

In  order  to  perform  all  the  necessary  multiply  operations,  we  secured  a  parallel  multiply  unit 
for  the  station  SDS  9300  computer.  This  unit  was  designed  to  perform  twelve  multiplications  in 
the  cycle  time  of  one.  Unfortunately,  delays  in  making  both  the  hardware  and  software  parts  of 
the  system  perform  properly  were  encountered,  and  a  cut  in  the  funds  available  to  support  the 
radar  site  forced  us  to  return  the  unit  to  the  manufacturer  before  it  could  be  brought  into  opera¬ 
tion  for  ionospheric  studies.  The  recent  development  of  the  Cooley-Tukey  algorithm  for  fast 
spectrum  analysis^  has  made  the  approach  described  here  considerably  more  practical,  and  a 
number  of  groups  have  now  constructed  fast  digital  processors  for  performing  spectrum  analysis 
calculations . 


3.  Second  Plan  for  New  System 

Recognizing  that  the  digital  autocorrelation  method  might  encounter  difficulties,  we  devised 
an  alternate  scheme  that  was  somewhat  less  ambitious  and  seemed  sure  to  work.  As  it  tran¬ 
spired,  this  second  plan  was  the  one  we  put  into  effect.  In  order  to  examine  several  height  in¬ 
tervals  simultaneously,  the  number  of  integrators  available  must  be  increased.  Therefore,  it 
seemed  logical  to  employ  the  digital  computer  to  perform  the  task  of  signal  integration,  while 
using  analog  elements  to  analyze  the  spectra  of  the  signals.  If  the  gated-receiver  approach  were 
used,  a  separate  filter  bank  would  be  required  for  each  height  examined.  However,  by  applying 
the  signals  continuously  to  a  bank  of  wide  filters,  acceptable  height  and  frequency  resolution  can 
be  achieved  when  the  filter  bandwidth  b  is  given  by 

b  =  I/T  Hz  (1) 

where  T  is  the  pulse  length  (seconds).  The  optimum  height  and  frequency  resolution  are  obtained 
when  the  power-vs -frequency  response  of  the  filter  is  matched  to  that  of  the  energy  in  the  pulse. 
That  is,  for  square  pulses  the  voltage  response  of  the  filter  should  be  of  the  form 


v(f)  -  v(fQ) 


sin  27rT(f  —  fQ) 
27rT(f  -  fQ) 


(2) 


where  f  is  the  filter  center  frequency.  When  this  is  the  case,  the  observed  spectrum  is  identi¬ 
cal  to  that  obtained  with  a  gated  receiver  and  very  narrow  filters.  Specifically,  the  measured 
power  at  a  frequency  f  is 


p(fQ) 


-r 


W(f)tj(f  -  fo) 


df 


(3) 
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where  W(f)t  is  the  fluctuation  spectrum  of  the  scattering  plasma  at  delay  t.  The  convolution 
function  j (f  —  f  )  is 


j(f-f  )  =  - £ - 

°  7rT(f  -  f  ) 


sin  27rT(f  —  fQ)i 
2irT(f-fo)  j 


(4) 


The  success  of  this  approach  depends  critically  upon  the  availability  of  filters  that  have 
characteristics  close  to  Eq.(2).  Fortunately,  multipole  crystal  filters  can  now  be  designed  using 
computer  programs,  so  that  suitable  filters  are  commercially  available. 

The  chief  difficulty  that  could  be  encountered  with  the  filter  bank  approach  is  that  it  is  neces¬ 
sary  to  take  a  large  number  of  samples  at  the  filter  outputs  simultaneously,  and  then  enter  these 
into  the  computer  serially.  The  problem  can  be  avoided  if  the  signals  applied  to  the  filters  are 
themselves  delayed  with  respect  to  each  other,  thereby  permitting  the  outputs  to  be  sampled  one 
at  a  time. 


4.  Realization  of  Spectrum  Analyzer 

Figure  1  provides  a  block  diagram  of  the  scheme  adopted  for  the  new  spectrum  analyzer. 

The  IF  signals  are  applied  to  a  tapped  delay  line.  Each  tap  drives  one  of  24  filters  that  form 
the  matched-filter  bank.  Arrangements  are  made  to  sample  the  filter  outputs  via  a  switch  (S^  in 
Fig.  1)  at  the  same  rate  as  the  signals  propagate  down  the  delay  line.  In  our  system,  an  acoustic 
delay  line  is  employed  and  the  pickup  coils  are  spaced  to  apply  replicas  of  the  IF  signals  to  the 
filters  with  a  relative  delay  of  20psec.  Switch  S^  consists  of  24  balanced-diode  gates  that  are 
opened  in  sequence,  one  at  a  time,  for  20|j.sec  by  a  clocked-digital  shift-register.  Thus,  in  or¬ 
der  to  sample  24  filters,  an  interval  of  480psec  is  required.  This  allows  a  separate  spectrum 
to  be  obtained  every  0.5  msec  along  the  timebase. 


200  kHl 
IF 

SIGNAL 


TAPPED  DELAY  LINE 


Fig.  1.  Block  diagram  of  new  spectrum  analyzer  scheme. 
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The  principal  drawback  of  this  scheme  is  that  a  separate  filter  bank  is  required  for  each  of 
the  pulse  lengths  used.  Thus,  two  new  filter  banks  were  purchased  to  suit  0.5-  and  1.0 -msec 
pulse  observations;  for  2.0 -msec  measurements,  a  third  bank  was  assembled  using  the  old 
500-Hz  filters  employed  in  the  original  spectrum  analyzer. 

5.  Filter  Networks 

Figure  2  shows  one  of  the  24  filter  networks  employed  in  the  spectrum  analyzer.  Each  pick¬ 
up  coil  on  the  delay  line  drives  a  pickup  amplifier  whose  gain  is  adjustable.  In  this  way,  the  at¬ 
tenuation  of  the  signals  down  the  delay  line  can  be  compensated.  The  pickup  amplifier  drives 
three  independent  amplifiers  that  are  coupled  to  the  three  filters  (matched  to  0.5-,  1.0-,  and 
2.0 -msec  pulses)  via  attenuators.  The  gains  of  these  amplifiers  differ  so  that  the  mean  noise 
voltage  appearing  at  the  filter  outputs  is  approximately  the  same  (i.e.,  independent  of  the  band¬ 
width).  The  filters  are  driven  continuously  and  a  relay-operated  switch  (S^  in  Fig.  2)  selects  the 
filter  in  use.  After  further  amplification,  the  signals  are  applied  to  the  diode  gates  that  act  as 
the  selector  switch  (S^  in  Fig.  1).  Each  diode  gate  is  normally  held  closed  by  a  flip-flop.  When 
allowed  to  open,  a  current  of  about  5  mA  flows  through  the  diodes.  This  limits  the  peak  RF  sig¬ 
nal  that  can  be  driven  through  the  gate  to  0.00  5R  volts,  where  R  is  the  output  resistance.  To 
reduce  the  effects  of  DC  offsets  and  switching  transients,  R  is  made  as  small  as  possible  (actu¬ 
ally  400  ohms)  consistent  with  being  able  to  obtain  10  volts  peak  detected  output.  The  diode 
bridges  can  be  balanced  to  yield  zero  output  when  no  signal  is  applied. 

6.  Frequency  Control 

In  order  to  obtain  suitable  filters,  we  were  obliged  to  employ  center  frequencies  in  the  vicin¬ 
ity  of  200  kHz.  Accordingly,  two  identical  banks  of  12  filters  were  obtained  for  each  pulse  length. 
These  filters  have  center  frequencies  spaced  at  2-kITz  intervals  and  span  the  interval  188  to 
210  kHz.  Each  filter  bank  is  driven  by  a  separate  receiver  channel,  as  shown  in  Fig.  3.  In  this 
way,  by  adjusting  the  2-MHz  to  200-kIlz  local-oscillator  (LO)  frequencies,  the  filters  can  be 
moved  to  cover  different  intervals.  For  UHF  observations,  the  filters  are  interlaced  to  yield 
measurements  at  l-kllz  intervals,  i.e.,  the  spectrum  analyzer  is  adjusted  to  examine  the  echoes 
at  — 11.5,  — 10.5,  —9.5,  .  .  .  ,—0.5,  +0.5,  .  .  .  ,  +11. 5 kHz  with  respect  to  the  radar  frequency.  The 
spectrum  analyzer  can  also  be  employed  with  the  E-band  radar  and,  in  this  case,  the  filter  banks 
are  placed  end  to  end,  i.e.,  the  range  —23  to  +23  kHz  with  respect  to  the  radar  frequency  is  ex¬ 
amined  at  2-kIIz  intervals.  Table  III  gives  the  LO  frequencies  required  for  these  settings. 

Figure  3  shows  that  the  odd-numbered  filters  (spanning  —11.5,  —9. 5,  .  .  .  ,  +10. 5 kllz  at  UHF) 
are  driven  with  replicas  of  the  signal  delayed  at  20,  40,  .  .  .  ,  240  psec.  In  order  that  the  even  fil¬ 
ters  (covering  — 10.5,— 8.5,  .  .  .  ,  +11. 5 kHz)  be  driven  at  later  delays,  a  fixed  length  of  delay  line 
introducing  240-psec  delay  is  employed.  The  attenuation  introduced  by  this  additional  piece  of 
line  is  35  dB,  requiring  that  an  additional  driver  amplifier  be  interposed  between  it  and  the  tapped 
delay  line  it  feeds. 

In  operation,  we  found  that  the  two  separate  mixers  (Fig.  3)  could  not  be  made  identical,  with 
the  result  that  the  spectra  obtained  yielded  points  that  were  alternately  high  and  low.  Despite 
several  attempts  at  improving  the  mixers  this  problem  could  not  be  eliminated,  and  we  decided 
to  commutate  the  mixer  outputs;  this  is  done  under  the  control  of  the  computer  every  m  sweeps 
of  the  timebase  (m  is  usually  set  at  750).  This  commutation  exchanges  the  roles  of  the  even  and 
odd  filters,  and  the  computer  keeps  track  of  the  correct  filter  center  frequencies  when  summing 
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Fig.  3.  Arrangement  of  delay-line  system. 
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TABLE  III 

LOCAL-OSCILLATOR  FREQUENCIES  EMPLOYED 

IN  THE  SPECTRUM  ANALYZER 

Mode 

Local 

Oscillator  1 
(MHz) 

Local 

Oscil lator  2 
(MHz) 

UHF 

Normal 

1 . 8005 

1.8015 

Transposed 

2.1995 

2.1985 

Normal 

1.789 

1.813 

L-band 

Transposed 

2.211 

2.187 

the  samples.  This  modification  successfully  eliminated  effects  due  to  dissimilarity  of  the  mixers. 
We  then  found  that  the  spectra  exhibited  an  asymmetry;  typically,  one  side  would  be  somewhat 
stronger  than  the  other.  This  is  attributable  to  unequal  responses  of  the  filters  in  their  skirts. 
Accordingly,  we  decided  to  transpose  the  filters  about  the  radar  frequency,  and  this  was  carried 
out  by  switching  the  LO  frequencies  from  below  2  MHz  to  an  equal  amount  above  (Table  III).  In 
effect,  the  filters  are  electrically  reflected  about  the  radar  frequency  since  they  are  employed 
to  examine  alternately  the  high  and  low  sidebands  of  the  spectrum.  Commutating  the  LO  frequen¬ 
cies  in  this  way  is  performed  every  2m  sweeps. 

These  changes  cause  each  filter  to  be  employed  at  four  different  places  in  the  spectrum; 
conversely,  each  point  in  the  spectrum  is  derived  from  samples  taken  from  four  physically  sep¬ 
arate  filters.  By  employing  the  computer  to  perform  the  commutating  and  keep  track  of  it,  the 
addition  of  the  samples  can  proceed  in  such  a  way  that  any  true  asymmetry  in  the  spectrum  is 
preserved.  That  is,  this  switching  serves  only  to  eliminate  instrumental  effects. 


7.  Rectifier  and  Sampling 

Figure  4  shows  the  rectifier  system  employed  to  sample  the  filter  outputs.  In  order  that  the 
rectifier  can  respond  within  20  psec  to  the  full  value  of  the  envelope  of  the  signal,  it  must  have  a 
wide  bandwidth  and  be  followed  with  a  time  constant  that  is  short  compared  with  20  psec.  Unfor¬ 
tunately,  as  noted  above,  we  were  obliged  to  purchase  filters  operating  at  ~200kHz  so  that  the 
rectifier  can  sample,  at  most,  four  cycles  of  the  waveform.  By  employing  a  full-wave  rectifier, 
the  period  of  the  ripple  applied  to  the  low-pass  filter  is  reduced  to  2.5 psec,  and  a  time  constant 
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Fig.  4.  Rectifier  system. 
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of  7  gsec  (the  geometric  mean  of  2.5  and  20  pscc)  was  found  to  be  an  acceptable  compromise.  As 
in  tli e  power  profile  channel,  a  linear  rectifier  is  employed  which  drives  a  10 -bit  analog-to- 
digital  converter  operating  over  a  10 -volt  range.  These  samples  of  signal  voltage  are  squared 
in  the  computer. 

We  found  that  the  electronic  switch  (Fig.  1)  introduced  noise  or  transients  whenever  the 
digital  shift-register  was  clocked.  These  transients  could  not  entirely  be  removed  by  filtering 
at  the  rectifier  output,  owing  to  the  requirement  that  the  time  constant  be  less  than  20  gsec. 
Considerable  improvement  was  obtained  by  filtering  at  the  input  (Fig.  4),  and  by  using  a  second 
diode  gate  to  connect  the  input  to  the  rectifier  to  ground  for  a  period  of  Bpsec  centered  on  the 
clock  time.  This  lowered  the  available  sample  time  to  15psec,  but  reduced  the  switching  noise 
to  the  level  of  the  least-significant  bit  in  the  A/D  converter.  The  5-psec  dead  period  also  is 
necessary  to  permit  the  low-pass  filter  to  discharge  between  samples,  otherwise  energy  from 
one  sample  can  be  carried  over  into  the  next.  Sampling  of  the  spectrum  channel  is  carried  out 
at  the  end  of  the  1  5-psec  sample  window,  i.e.,  immediately  before  the  input  to  the  rectifier  is 
shorted  (by  in  Fig.  4).  The  timing  of  this  chain  of  sample  pulses  is  adjusted  to  allow  for  the 
overall  delay  in  the  receiver,  filter,  and  rectifier  to  yield  samples  that  are  precisely  controlled 
in  delay  with  respect  to  the  leading  edge  of  the  transmitter  pulse.  That  is,  to  obtain  a  spectrum 
at  a  nominal  delay  of  1.5  msec,  the  first  sample  is  taken  approximately  2.1  msec  following  the 
transmitter  pulse. 

G.  Timing  Equipment 

All  events  in  the  cycle  of  operations  are  timed  by  digital  clocks  operating  from  a  100 -kHz 
standard  frequency  derived  from  the  station  master  oscillator.  Control  of  the  transmitter  pulse 
length,  sweep  interval,  receiver  gate  position,  location,  and  duration  of  the  noise  calibration 
pulse  are  under  control  of  a  radar  timer  that  provides  four  preset  options  or  modes.  These 
modes  are  designated  A  through  I)  and  are  selectable  from  the  control  console  by  push  buttons. 
Alternatively,  the  timer  can  be  placed  under  remote  control,  in  which  case  the  mode  selection 
is  performed  by  the  computer. 

III.  OPERATING  PROCEDURE 
A.  Operating  Modes 

As  outlined  in  the  previous  section,  the  new  equipment  provides  four  selectable  operating 
modes  designated  A  through  D.  Relay  controls  select  the  appropriate  bandwidth  filters  for  each 
mode,  and  the  time  constant  of  the  low-pass  filter  in  the  power  profile  channel.  Table  IV  sum¬ 
marizes  the  operating  parameters  for  each  of  these  modes,  and  the  location  on  the  timebase  of 
the  noise  calibration  pulse  and  its  width.  Also  given  is  the  position  and  width  of  a  portion  of  the 
timebase,  beyond  all  ionospheric  echoes  and  preceding  the  noise  calibration  pulse,  that  is  sam¬ 
pled  as  a  measure  of  the  intrinsic  receiver  noise.  In  the  spectrum  channel,  a  separate  spectrum 
is  available  every  0.5msec.  These  are  completely  independent  only  when  0.5-msec  pulses  are 
employed  (mode  B),  but,  despite  this,  all  possible  spectra  are  accepted  between  delays  on  the 
timebase  given  in  Table  IV.  An  additional  five  independent  spectra  (i.e.,  spaced  by  the  pulse 
length)  are  taken  at  large  delays  on  the  timebase  (see  table)  to  measure  the  noise  present  in 
each  channel. 
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PARAMETERS  EMPLOYED 

TABLE  IV 

FOR  THE  FOUR  SELECTABLE  MODES 

Mode  (msec) 

Porometer 

A 

B 

C 

D 

Pulse  length 

0.1 

0.5 

1.0 

2.0 

Interpulse  period 

8.0 

13.0 

20.0 

40.0 

Receiver  gate  deloy* 

0.5 

1.5 

2.0 

5.0 

Profile  chonnel  sample  spacing 

0.05 

0.20 

0.20 

0.20 

Profile  channel  time  constont 

0.05 

0.20 

0.20 

0.20 

Deloy  ta  receiver  naise  window 

5.5 

9.0 

12.0 

17.0 

Width  af  receiver  naise  window 

0.5 

1.0 

2.0 

2.0 

Delay  ta  calibration  naise  pulse 

6.5 

11.0 

17.0 

35.0 

Width  af  calibrotian  naise  pulse 

0.5 

1.0 

2.0 

2.0 

Deloy  ta  first  signal  spectrum 

- 

1.5 

3.0 

6.0 

Delay  ta  lost  signal  spectrum 

- 

5.0 

7.5 

10.5 

Delay  ta  first  naise  spectrum 

- 

8.0 

8.0t 

11.0+ 

Delay  ta  lost  naise  spectrum 

- 

10.5 

13.0t 

21.0+ 

*  May  be  manually  adjusted. 

t  Increosed  by  2.0  msec  in  February  1970. 

B.  Real-Time  Computer  Operations 

A  computer  program  has  been  developed  that  affords  complete  control  over  the  experiment. 
Basically,  the  program  squares  and  sums  the  signal  samples  available  at  the  outputs  of  the  power 
profile  and  spectrum  channels.  In  the  case  of  the  power  profile,  it  is  necessary  only  to  keep 
track  of  the  positions  of  the  samples  along  the  timebase.  However,  the  spectrum  channel  pro¬ 
vides  four  separate  sample  sequences,  depending  upon  whether  or  not  the  mixer  outputs  have 
been  interchanged  and  the  LO  frequencies  are  above  or  below  the  sample  frequency.  For  this 
reason,  the  computer  controls  these  switching  operations  via  an  output  register  and  associated 
relays.  The  interval  between  commutations  is  determined  by  counting  a  predetermined  number 
of  sweeps,  and  the  duration  of  the  entire  integration  period  is  determined  in  the  same  fashion. 

All  the  parameters  listed  in  Table  IV,  as  well  as  the  duration  of  the  integration  period,  am¬ 
plitude  of  the  calibration  pulse,  commutation  interval,  date,  etc.,  are  stored  in  the  program  in 
a  manner  that  makes  them  accessible  to  change  via  the  typewriter. 

When  the  radar  timer  is  placed  in  remote  control,  the  computer  takes  over  the  selection  of 
the  operating  modes  via  an  output  register.  Two  options  are  then  available.  In  one,  the  opera¬ 
tor  requests  a  run  of  a  given  type  (A  to  D)  by  pressing  the  appropriate  key  on  the  typewriter;  the 
computer  then  sets  the  radar  into  that  mode  and  commences  taking  data  at  the  next  round  10 -sec 
interval.  The  alternative  is  for  the  operator  to  elect  an  "automatic  sequence"  option  at  the  com¬ 
mencement  of  the  operating  period.  In  this  case,  it  is  necessary  to  specify  the  cycle  in  which 
runs  are  to  be  carried  out  ( i . e . ,  the  type  and  duration  of  each  run  before  the  cycle  repeats). 
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Automatic  sequencing  may  be  interrupted  to  reject  a  particular  run,  change  the  stored  param¬ 
eters,  or  change  the  desired  sequence.  At  the  commencement  of  the  first  run  taken  following 
Greenwich  midnight,  the  computer  requests  the  operator  to  supply  a  new  date. 

C.  F-Region  Critical  Frequency 

An  important  parameter  in  the  data  processing  (Sec.  IV-C)  is  the  critical  frequency  of  the 
F-region  determined  using  a  C-4  ionosonde  located  about  1km  from  the  radar.  In  the  past,  this 
sounder  has  been  left  to  operate  unattended,  and  records  were  made  on  35-mm  film.  In  many 
instances,  during  radar  operations,  malfunctions  occurred  that  were  not  discovered  at  the  time. 

In  order  to  remedy  this,  arrangements  have  been  made  to  permit  the  ionosonde  to  be  turned  on 
from  the  radar  console,  and  control  of  its  operating  frequency  transferred  to  a  local  frequency 
synthesizer.  The  ionosonde  receiver  output  may  also  be  viewed  on  a  local  A-scope  display.  In 
this  way,  the  operator  may  measure  fQF2  remotely  by  advancing  the  frequency  of  the  synthesizer 
until  the  delay  to  the  ordinary  trace  increases  rapidly.  These  measurements  of  the  critical  fre¬ 
quency  are  stored  in  the  computer,  which  requests  the  value  of  fQF2  each  time  an  A-mode  run 
is  reached  in  the  mode  sequence,  or  is  requested  by  the  operator.  Film  records  of  automatic 
sweeps  made  by  the  ionosonde  at  30 -minute  intervals  are  also  made,  and,  if  desired,  values  of 
f^F2  obtained  from  these  may  be  used  in  the  data  analysis  in  place  of  those  obtained  by  the  radar 
operator. 

D.  Computer  Output 

At  the  end  of  the  run  period,  the  computer  halts  the  integration  process  and  partially  proc¬ 
esses  the  results  before  proceeding  to  the  next  run.  These  data  are  written  on  magnetic  tape 
for  later  non-real -time  analysis  (Sec.  IV),  and  are  printed  by  a  high-speed  line  printer  to  allow 
the  operator  to  see  that  the  experiment  is  proceeding  satisfactorily.  In  the  case  of  the  power 
profile  results,  the  raw  computer  sums  are  tabulated  vs  delay  t  along  the  timebase  and  presented 
on  a  monitor  display.  The  equivalent  ionospheric  height  h  is  computed  according  to 

h  =  150(t  -  km  (5) 

where  T  is  the  transmitter  pulse  length,  and  t  is  the  total  receiver  delay  (all  times  being  in 
milliseconds).  The  computer  also  calculates  the  mean  of  the  samples  within  the  receiver  noise 
window  and  subtracts  this  from  all  other  samples.  The  residuals  are  then  converted  to  echo  tem¬ 
perature  T^  employing  the  known  amplitude  (°K)  of  the  calibration  pulse  as  scale  factor,  Finally, 
an  echo  power  profile  is  computed  according  to 

a  =  2T  h2kb  W  km2  (6) 

r  r 

where  k  is  Boltzmann's  constant  (1.38  X  10”  erg/°K),  and  b  is  the  receiver  predetector  band¬ 
width  (in  hertz).  The  quantities  h,  Tr#  and  c are  also  tabulated  against  delay  t  on  the  high¬ 
speed  printout. 

It  may  be  shown  that  cr  is  directly  related  to  the  echo  cross  section  per  unit  height  in  the 

r  i 

ionosphere,  and  hence  proportional  to  the  electron  density.  Accordingly,  at  the  completion  of 
each  run,  the  computer  also  presents  on  the  monitor  display  a  plot  of  o ^  normalized  to  the  larg¬ 
est  value. 
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For  the  spectrum  channel,  the  computer  calculates  the  mean  noise  intensity  Pn(f)  in  each 
channel  and  subtracts  this  from  the  sums  at  the  other  delays.  Finally,  these  residuals  are  di¬ 
vided  by  Pn(f)  to  yield  the  signal-to-noise  ratio  at  each  delay  and  frequency.  These  values  are 
printed  out  and  stored  on  the  magnetic  tape. 

E.  Observing  Program 

The  new  equipment  provides  the  opportunity  for  making  observations  in  a  variety  of  different 
ways.  However,  since  each  possible  observing  procedure  may  require  a  separate  program  to 
handle  the  data  analysis,  we  have  thus  far  limited  the  observations  to  two  kinds.  The  first  of 
these  provides  good  coverage  of  the  entire  altitude  interval  150  to  1100  km  with  modest  time 
resolution  (30  minutes).  For  this  program,  A-,  B-,  and  C-mode  runs  are  performed  sequen¬ 
tially,  each  being  8  minutes  in  duration.  A  second  type  of  observing  schedule  has  been  employed 
to  emphasize  the  drift  measuring  capability  afforded  by  the  C-mode  runs.  For  these  observa¬ 
tions,  the  sequence  of  runs  is  A,  B  (each  4  minutes),  followed  by  four  C-mode  runs  (each  lasting 
8  minutes).  This  cycle  is  termed  "Drifts"  to  distinguish  it  from  the  normal  or  "Regular"  cycle. 

We  have  not  found  it  profitable  to  include  D-mode  runs  in  either  of  these  sequences,  owing 

to  the  low  signal-to-noise  ratio  encountered  at  the  altitudes  covered  with  this  mode  (Table  IV). 

However,  separate  observations  using  only  the  D-mode  are  planned  to  search  for  the  transition 

-f  -f 

altitude  between  O  and  H  ions,  which  has  not  been  detected  in  the  Millstone  observations  car¬ 
ried  out  thus  far. 

During  the  latter  part  of  1968,  we  attempted  to  make  observations  with  both  the  Regular 
and  Drifts  cycles  twice  per  month  for  periods  of  24  hours.  In  1969,  this  was  reduced  to  one 
24-hour  operation  of  each  type  per  month.  Periods  when  these  measurements  were  carried  out 
are  listed  in  Tables  V  and  VI. 

IV.  ANALYSIS 

A.  Introduction 

A  computer  program,  which  employs  as  input  the  digital  data  tape  on  which  all  the  radar 
data  are  stored,  has  been  written  to  analyze  the  results  obtained  in  the  Regular  and  Drifts  obser¬ 
vations.  Early  in  the  development  of  the  new  system,  we  planned  to  perform  all  data  processing 
at  the  completion  of  each  run,  or  cycle  of  runs,  thereby  eliminating  any  post-processing  and 
providing  the  final  results  in  usable  form  in  real  time.  This  still  appears  to  be  a  desirable  goal, 
but  because  the  data  collection  and  data  analysis  programs  both  grew  in  complexity,  we  decided 
to  separate  them,  allowing  data  collection  to  commence  prior  to  the  availability  of  an  entirely 
satisfactory  analysis  program. 

The  analysis  program  described  here  has  been  under  development  for  some  time  and,  al¬ 
though  open  to  further  refinement,  is  now  probably  close  to  final  form.  The  program  provides 
a  smooth  electron  density  profile  that  is  corrected  for  the  effects  of  the  ratio  of  T^/T.  on  the 
electron  scattering  cross  section,  and  yields  electron  temperatures  T^  corrected  for  the  effects 
of  the  variation  of  Debye  length  with  altitude.  As  such,  the  density  measurements  depend  upon 
the  temperature  measurements,  and  vice  versa,  so  that  an  iterative  processing  scheme  is 
required. 

Table  VII  lists  the  data  provided  by  the  programs  discussed  in  the  previous  section.  As  is 
evident  from  the  table,  the  height  resolution  for  the  density  measurements  considerably  exceeds 
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TABLE  V 

REGULAR  UHF  OBSERVATIONS 

PRIOR  TO  JANUARY  1970 

Dote 

Start  (GMT) 

End  (GMT) 

Comments 

1968 

15-16  August 

1500 

1400 

22-23  August 

1800 

1700 

29-30  August 

1300 

1300 

12-13  September 

2000 

2000 

19-20  September 

1300 

1300 

3-4  October 

2030 

2230 

22-23  October 

1300 

1300 

14-15  November 

1530 

1530 

25-26  November 

1500 

1500 

17-18  December 

1400 

1400 

30-31  December 

1500 

1530 

1969 

16-17  Jonuory 

1430 

1430 

5-6  February 

1430 

1400 

26-27  February 

1530 

1530 

Antenna  gain  low  owing 
to  snow 

21-22  Morch 

2120 

1010 

23-24  March 

0140 

1000 

Auroro 

9-10  April 

1430 

1430 

6-7  May 

1300 

1300 

5-6  June 

1330 

1400 

1-2  July 

1300 

1315 

14-15  August 

1315 

1345 

9-10  September 

1415 

1400 

30  September,  1  October 

0300 

1400 

Auroro 

3-4  November 

1900 

2200 

8-9  December 

1745 

1745 
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TABLE  VI 

DRIFTS 

UHF  OBSERVATIONS  PRIOR  TO  JANUARY  1970 

Date 

Start  (GMT) 

End  (GMT) 

Comments 

1968 

2-3  August 

0100 

0200 

Na  cammutatian 

8-9  August 

0230 

1130 

Na  commutation 

1 3-14  August 

1530 

1530 

Na  cammutatian 

20-21  August 

1300 

1400 

27-28  August 

1330 

1330 

2-3  September 

1700 

1700 

11-12  September 

1800 

1800 

26-27  September 

1400 

1400 

8-9  October 

1600 

1800 

29-30  October 

1500 

1500 

1  9  Navember 

1700 

2300 

29-30  Navember 

1400 

1400 

12-13  December 

1700 

1800 

1969 

31  January 

0100 

1300 

12-13  February 

1420 

1450 

25-26  March 

2130 

2130 

23-24  April 

0320 

0310 

30  May 

1300 

2200 

Commutating  failed 
at  2200 

1  June 

0050 

1200 

Bad  data  tape  prior 
ta  1  June 

23-24  June 

1320 

1310 

9-10  July 

1140 

0945 

29-30  July 

2300 

2400 

26-27  August 

2030 

2030 

23-24  September 

1300 

1300 

14-15  October 

2200 

2200 

20-21  Navember 

1330 

1500 

22-23  December 

2045 

2045 
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TABLE  VII 

data  AVAILABLE  FOR  PROCESSING 

Height 

Resolution 

Sample 

Spacing 

(Density) 

Heights  for  Temperature 
Measurements 

Mode 

(km) 

(km) 

(km) 

A 

15 

7.5 

None 

B 

75 

30 

225,  300,  375,  450,  525, 

600,  675 

C 

150 

30 

450  ,  525  ,  600  ,  675  ,  750, 

825,  900,  975,  1050,  1125 

that  for  the  temperature  measurements.  Advantage  is  taken  of  this  when  constructing  the  den¬ 
sity  curve  by  using  only  the  A-mode  data  up  to  an  altitude  above  hmaxF2. 

B.  Electron  and  Ion  Temperatures 

In  essence,  the  electron  and  ion  temperatures  are  obtained  by  comparing  the  observed  spec¬ 
tra  with  those  computed  theoretically  for  different  temperatures.  These  theoretical  spectra  were 
calculated  using  expressions  published  by  Fejer,  and  assume  the  following:  the  electrons  and 
ions  each  have  a  Maxwellian  velocity  distribution,  but  need  not  be  at  the  same  temperature;  all 
ions  are  singly  charged;  collisions  with  neutrals  and  the  influence  of  the  earth's  magnetic  field 
may  be  neglected.  For  the  observations  described  here,  we  believe  all  these  conditions  have 

A  A  -f 

been  met.  In  carrying  out  the  calculations,  we  also  assume  that  O  is  the  only  ion  present. 

This  assumption  appears  to  be  valid  at  Millstone,  at  least  over  most  of  the  altitudes  examined. 

For  the  measurements  made  at  225-km  altitude,  the  pulse  also  illuminates  regions  beneath  where 
and  NO  are  present  —  especially  in  summer.  '  Neglecting  the  presence  of  these  ions  in 

the  reduction  should  cause  both  T  and  T  to  be  slightly  underestimated. 

e  1  +  + 

At  high  altitudes  (>7  50  km),  it  is  possible  that  small  numbers  of  light  ions  (II  ,  He  )  are  pres¬ 
ent  which  will  cause  T.  to  be  overestimated.  However,  from  the  observed  signal  spectra,  we 

*  4. 

have  not  yet  been  able  to  identify  regions  in  which  a  significant  abundance  (^25  percent)  of  H  or 
He^  is  present.  A  final  assumption  made  in  the  calculations  is  that  the  plasma  frequency  is 
10  MHz.  This  is  equivalent  to  assuming  that  the  Debye  length  is  very  small  and  is  justified 
only  at  altitudes  below  about  450  km.  Once  an  estimate  has  been  made  of  the  Debye  length,  it  is 
possible  to  correct  for  this  assumption  as  discussed  below. 

The  fluctuation  spectrum  W(f)  of  the  plasma,  calculated  in  the  manner  outlined  above,  is 
next  convolved  with  the  function  [Eq.  (4))  representing  the  instrumental  resolution  set  by  the  pulse 
length  T.  With  the  computer  program  presently  available,  these  calculations  do  not  proceed 
quickly  enough  to  permit  a  point -by -point  comparison  of  the  measured  and  theoretical  spectra 
in  a  least-mean-squares  sense.  Instead,  the  reduction  is  carried  out  by  comparing  two  param¬ 
eters  that  characterize  the  spectra  —  namely,  the  ratio  of  the  power  in  the  wings  to  that  at  the 
center  frequency,  and  the  spectrum  width  (to  a  point  of  half-peak  power).  The  ratio  and  width 
of  the  theoretical  spectra  were  obtained,  and  polynomial  expressions  were  then  derived  (by  least- 
mean-squares  fitting  in  a  computer)  to  yield  Te  and  T.  as  functions  of  these  two  parameters. 
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These  functions  each  employ  ten  terms  and  allow  the  measured  ratio  and  width  to  be  converted 
to  temperature  to  an  accuracy  of  1  percent  or  better.  Different  functions  are  required  for  each 
of  the  pulse  lengths  employed,  and,  in  order  to  obtain  sufficient  accuracy,  different  equations 
are  employed  depending  upon  whether  Tg/T.  is  greater  or  less  than  2.0. 

To  recover  the  ratio  parameter  for  each  measured  spectrum,  parabolas  are  fitted  to  the 
four  center  points  and  four  points  in  the  wing  of  the  spectrum  chosen  to  include  the  highest  adja¬ 
cent  pair.  Figures  5  and  6  are  flow  diagrams  showing  the  procedure  employed  for  fitting  the 
two  parabolas  and  ensuring  that  they  are  of  correct  shape.  In  order  to  locate  the  position  in  fre¬ 
quency  of  the  half-peak  value,  linear  interpolation  is  employed  between  the  pair  of  points  on  the 
edge  of  the  spectrum  whose  amplitudes  bracket  the  half-peak  power  value. 


Fig.  5.  Method  of  renormalizing  by  fitting 
a  parobolo  to  center  of  spectrum. 


When  determining  the  ratio  and  width  parameters,  the  two  halves  of  the  spectrum  are  an¬ 
alyzed  separately,  and  independently  converted  to  estimates  of  electron  and  ion  temperatures. 

In  addition,  a  mean  is  taken  of  the  ratios  and  widths  and  is  also  used  to  obtain  T  and  T..  Be- 

e  i 

cause  it  should  be  insensitive  to  shifts  in  the  spectrum  introduced  by  drifts  (Sec.  IV-D),  this  mean 
value  is  adopted  as  the  best  estimate  where  both  halves  of  the  spectrum  have  been  successfully 
analyzed. 

Moorcroft*4  has  shown  that  the  effect  of  increasing  the  Debye  length  D^  on  the  spectrum 
shape  is  exactly  equivalent  to  lowering  the  ratio  Te/T..  Thus,  spectra  interpreted  on  the  assump¬ 
tion  that  D^  —  0  yield  a  fictitious  value  /?  for  the  electron -to-ion-temperature  ratio  that  is  related 

to  the  true  value  T  /T.  in*4 
e'  l 


^  =  (1  +  a2)  /?  1  (7) 

i 

2  l/2 

where  a  -  47rDe/X,  and  Dq  =  (kTe/47rNe  j  '  c  in  which  N  is  the  electron  density  and  e  is  the 
electronic  charge.  Since  the  estimate  of  T.  does  not  depend  upon  o,  this  means  that  only  the 
initial  estimate  T^  of  the  electron  temperature  is  in  error.*4  For  our  wavelength,  the  correction 
formula  reduces  to 


(<-^r 


T  = 
e 


(8) 
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Method  of  finding  height  of  wing  by  fitting  a  parabola. 


Evidently,  to  obtain  the  true  temperature,  it  is  necessary  to  know  the  density  N.  Thus,  at 
this  stage,  reduction  of  the  temperature  data  is  halted  until  the  density  profile  has  been  reduced. 
When  this  has  been  done,  corrected  values  of  are  obtained  via  Eq.  (8).  These  are  next  em¬ 
ployed  in  a  recomputation  of  the  density  profile,  and  this  revised  estimate  of  N  is  then  used  to 
yield  a  second,  more  accurate  estimate  of  T  .  We  feel  that  the  accuracy  achieved  this  way  is 
sufficient  for  all  practical  purposes,  but,  of  course,  further  iterations  could  be  performed. 

C.  Electron  Density 

The  electron  density-vs-altitude  distribution  is  obtained  by  combining  data  obtained  in  the 
A-,  B-,  and  C-mode  runs  to  yield  the  shape  of  the  profile,  and  by  normalizing  this  to  have  the 
correct  peak  density  set  by 

N  =  1.24  x  10^(f  F2)2  electrons/cm^  (9) 

max  o  ' 

where  fQF2  is  in  megahertz.  The  value  of  the  F-region  critical  frequency  is  measured  and  stored 
prior  to  each  A-mode  run  (Sec.  III-C),  hence,  all  the  information  needed  to  process  the  results 
is  available  on  the  data  tape.  However,  a  new  table  of  f^F2  values  may  be  substituted  via  punched 
cards. 

The  principal  difficulty  encountered  in  writing  the  present  computer  program  was  in  devising 
a  way  of  combining  the  A,  B,  and  C  results  to  obtain  the  shape  of  the  profile.  Previously,  the 
parameter  log^a^  (Eq.  (6)]  was  machine  plotted,  and  a  single  curve  was  obtained  by  overlaying 
these  curves  for  the  best  fit.  That  is,  the  A-mode  results,  which  provide  the  best  height  resolu¬ 
tion,  were  allowed  to  control  the  results  up  to  some  point  above  hmax  where  they  became  noisy. 

A  smooth  transition  was  then  made  to  the  B-mode  results,  and  at  some  higher  altitude  to  the 
C  data.  This  combined  curve  was  then  corrected  for  the  effect  on  the  scattering  cross  section 

of  T  /T.  ( lief.  1). 

e'  i  ' 

When  attempting  to  obtain  a  computer  program  to  perform  these  tasks,  we  encountered  con¬ 
siderable  difficulty  in  devising  a  method  that  yielded  a  smooth  profile.  Typically,  changes  in 
the  density  gradient  would  be  encountered  whenever  there  was  a  transition  from  the  data  obtained 
in  one  mode  to  another.  In  the  method  finally  adopted,  the  data  are  renormalized  so  that  in  cer¬ 
tain  altitude  ranges  common  to  two  modes  the  points  have  a  common  center-  of  gravity.  A  single 
curve  through  all  the  points  is  then  obtained  by  least-mean-squares,  each  point  being  weighted 
in  proportion  to  its  signal-to-noise  ratio. 

The  program  proceeds  by  treating  all  the  data  obtained  between  successive  A-mode  runs  as 
a  set.  Where  more  than  one  run  (B,  C,  or  D)  of  the  same  kind  has  been  made  (as  in  the  Drifts 
program)  the  values  of  the  echo  power  T^  are  summed  and  divided  by  \fnq ,  where  n  is  the  num¬ 
ber  of  runs  of  the  same  type,  and  q  is  the  ratio  of  their  length  to  the  length  of  the  A-mode  run. 
This  serves  to  normalize  the  results  so  that  the  relative  uncertainty  in  each  point  is  inversely 
proportional  to  its  absolute  value. 

Owing  to  the  fact  that  the  density  falls  by  about  two  orders  of  magnitude  in  the  region  between 
hmax  anc*  ^00  km,  it  becomes  difficult  to  employ  a  least -mean-squares  curve-fitting  procedure, 
since  the  shape  of  the  curve  will  be  dominated  by  the  large  values  of  density  at  low  altitudes. 

For  this  reason  and  because  the  density  decreases  exponentially,  or  very  nearly  so,  it  is  con¬ 
venient  to  work  with  the  logarithm  of  the  density.  Thus,  the  logarithm  is  taken  of  all  the  values 
of  Tr  obtained  in  each  run  over  the  height  range  for  which  Tr  >  0.  Next,  the  effects  of  height  h 
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and  temperature  ratio  on  the  results  are  removed  by  computing 

logl0N  =  g(h)  +  const  (10) 

where 

g(h)  =  log10  (Tr)  +  2  log1()  (h)  +  log1Q  (1  +  P)  (11) 

and  p  is  the  initial  estimate  of  Tg/Tj  [Eq.  (7)]. 

To  obtain  the  height  dependence  of  /3  required  in  Eq.  (11),  the  following  scheme  is  employed. 
For  h<  22  5km,  linear  interpolation  between  the  value  observed  at  22  5km  and  an  assumed  value 
/ 3  =  1.0  at  130  km  (Ref.  12)  is  employed.  Between  22  5  and  600  km,  all  the  available  values  of  /3 
at  each  altitude  (Table  VII)  are  averaged  and  /3(h)  again  obtained  by  linear  interpolation  between 
each  pair  of  points.  A  straight  line  is  fitted  by  least -mean-squares  to  all  the  values  obtained  in 
the  C-mode  for  h  ^  600  km  and  employed  at  these  altitudes  for  the  height  dependence  of  /3,  pro¬ 
vided  that  /3  —  1.0  as  h  —•  °°.  In  the  event  that  the  slope  of  the  fitted  straight  line  is  positive,  lin¬ 
ear  interpolation  is  used  between  the  observed  mean  value  at  600  km  and  an  assumed  value  /3  =  1.0 
at  1500  km.  These  procedures  serve  to  ensure  that  the  third  term  in  Eq.  (11)  is  well  behaved. 

The  discontinuities  in  the  gradient  d/3/dh  that  are  introduced  cause  only  minor  changes  in  the 
density  curve;  these  are  largely  removed  by  the  subsequent  smoothing. 

The  next  step  is  to  renormalize  the  g(h)  values  for  the  A,  B,  and  C  data  so  that,  for  certain 
common  altitude  regions  above  hmax>  they  have  the  same  centers  of  gravity.  To  do  this,  use  is 
made  of  the  fact  that  dg(h)/dh  is  nearly  constant.  Thus,  a  curve  of  second  order  is  fitted  to  the 

g(h)  results  for  the  A-mode  data  above  h  +  60  km.  A  similar  curve  is  fitted  to  the  B-mode 
6  max 

data  above  h  +  200km.  A  shift  is  then  made  to  the  B-mode  curve  to  cause  it  to  intersect  the 
max 

A-mode  curve  at  the  midpoint  of  the  height  interval  common  to  the  two  sets  of  data.  The  C-mode 
results  are  normalized  to  the  B-mode  results  in  a  similar  fashion. 

In  fitting  these  second-order  functions  to  the  data,  an  attempt  is  made  to  remove  spurious 
points  caused,  for  example,  by  echoes  from  satellites  that  occasionally  pass  through  the  beam. 
Thus,  when  the  curve  is  first  fitted  each  point  is  weighted  inversely  as  the  square  of  its  height. 
Next,  all  points  that  deviate  from  the  curve  by  more  than  five  standard  deviations  are  removed, 
and  the  curve  refitted  to  the  remaining  points.  In  the  second  fit,  each  point  is  weighted  by  the 
original  value  of  temperature  T  .  (Although  the  relative  uncertainty  in  the  points  is  proportional 

_  i  r 

to  T^  ,  the  temperature  cannot  be  used  to  weight  the  points  until  spuriously  large  ones  have  been 
removed. ) 

With  the  three  sets  of  points  thus  normalized,  a  new  second-order  function  is  fitted  to  all 
the  points  between  h max  +  60  km  and  the  highest  point  available,  each  point  being  weighted  ac¬ 
cording  to  its  value.  This  curve  is  taken  as  the  variation  of  g(h)  above  h  +  100  km,  and  below 
this  altitude  the  unsmoothed  A-mode  data  are  employed.  The  constant  term  in  Eq.  (10)  is  deter¬ 
mined  by  fitting  a  parabola  to  the  points  in  the  vicinity  of  the  peak  and  renormalizing  so  that 

N  =  N  [Eq.  (9)]  at  h 

max  1  M  1  max 

These  steps  yield  an  initial  estimate  of  the  density  variation  that  is  employed  to  correct  the 
electron  temperature  results  via  Eq.{8).  As  discussed  above,  this  yields  revised  values  of  the 
ratio  T^/T.  permitting  a  new  estimate  of  g(h)  from 

g(h)  =  log10Tr  +  2  1og10h+log10[l  +  (Te/T.)J  .  (12) 
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The  variation  of  Te/T.  with  altitude  is  obtained  in  exactly  the  same  way  as  the  variation  of  the 
initial  estimate  /?.  The  curve  fitting,  renormalization,  and  final  smoothing  are  then  all  per¬ 
formed  again,  exactly  in  the  same  manner,  to  yield  the  final  density  curve. 


D.  Vertical  Drift  Velocity 

If  the  ionospheric  plasma  in  the  scattering  volume  sampled  by  the  radar  has  a  mean  radial 
velocity  v^  with  respect  to  the  radar,  the  entire  spectrum  will  be  shifted  by  an  amount 


Afd  = 


~2vd 

X 


Hz 


(13) 


For  the  Millstone  Hill  radar  (X  =  68  cm),  this  corresponds  to 

Af^ - 3v^  Hz  (14) 

when  v^  is  in  meters/second. 

Since  the  vertical  drifts  in  the  ionosphere  are  at  most  a  few  tens  of  meters  per  second,  Af^ 
is  typically  less  than  100  Hz.  The  overall  width  of  the  spectrum  on  the  other  hand  usually  exceeds 
10  kHz,  thus,  extremely  precise  spectrum  measurements  are  required  if  the  presence  of  a  drift 
of  the  plasma  is  to  be  detected.  Despite  this,  the  results  of  the  C-mode  observations  clearly 

reveal  effects  due  to  drifts.  Various  methods  of  extracting  the  drift  information  have  been  tried 

i  s 

and  compared,  and  the  results  have  been  published  elsewhere.  Here,  we  shall  simply  sum¬ 
marize  the  approach  that  was  finally  taken. 

Woodman  and  Hagfors^  have  discussed  in  detail  the  method  and  accuracy  of  the  velocity  de¬ 
termination  procedure  adopted  at  Jicamarca.  In  this  scheme,  the  average  phase  shift  A cp  between 
a  pair  of  short  transmitted  pulses  separated  by  r  seconds  is  obtained.  Since  the  phase  shift  A (p 
is  determined  by 

A (p  -  27rAf^r  radians  (15) 


it  is  possible  to  recover  the  drift  velocity.  The  accuracy  of  measurement  initially  improves  with 
increasing  the  pulse  separation  r,  but  subsequently  falls  because  the  pulses  become  uncorrelated. 
Thus,  Woodman  and  Hagfors^  give  for  the  rms  velocity  uncertainty 


Avd- 


1  _L  f  1  -S(t)211/2 
(2n)l/2  4v  tS(t) 


m/sec 


(16) 


where  n  is  the  number  of  pulse  pairs  transmitted,  and  S(t)  is  the  normalized  echo  amplitude 
correlation  for  pulses  spaced  t  sec  apart. 

The  great  advantage  of  the  complex  correlation  function  approach  is  that  (for  symmetrical 
spectra)  the  phase  of  the  function  contains  all  (and  only)  the  information  concerning  the  doppler 
shift,  while  the  shape  is  defined  separately  by  the  normalized  echo  amplitude  correlation  func¬ 
tion  S(t).  These  advantages  can  be  achieved  for  the  Millstone  Hill  data  by  taking  the  complex 
Fourier  transform  of  the  measured  power  spectrum.  A  set  of  real  a(T>  and  imaginary  b(r)  parts 
are  then  obtained  where 


r  =  0,  1,  .  .  .  ,  ll(NF)’*  sec 

in  which  N  is  the  number  of  points  in  the  original  spectrum  (24),  and  F  is  their  spacing  (1kHz). 
The  phase  shift  A at  each  of  these  spacings  is  simply 
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(17) 


In  order  to  use  all  the  possible  estimates  of  A (p^,  two  different  approaches  have  been  tried.  In 
the  first,  a  least-mean-squares  straight  line  was  fitted  to  the  A (p ^  results.  However,  only  the 
first  6  of  the  12  points  were  considered  significant,  and  each  of  these  was  first  weighted  accord¬ 
ing  to  the  square  of  the  amplitude  correlation  at  that  spacing,  i.e.,  an  amount 


c2(r)  a2(r)+b2(r) 


(18) 


The  slope  a  of  this  least-mean-squares  straight  line  was  taken  to  be  the  rate  of  change  of  phase 
with  spacing  A<p  /r  and  employed  in  Eq.  (15)  to  determine  v^. 

The  least-mean-squares  approach  was  later  recognized  to  be  "suboptimum ,"  since  points  at 
large  spacings  r  tend  to  control  the  slope  cv.  It  may  be  shown  that  the  optimum  procedure  is 


to  search  for  the  value  of  o'  that  maximizes  the  sum^ 


(19) 


T 


in  which  a(T)  and  b(T)  are  the  real  and  imaginary  parts  of  the  transformed  measured  power  spec¬ 
trum,  S(t)  is  the  expected  correlation  (root  sum  square  of  real  and  imaginary  parts)  at  each 
spacing,  and  a  is  the  unknown  slope  of  the  phase  function,  i.e.,  A (p/r.  This  process  amounts 
to  multiplying  each  vector  (a  +  jb)  by  a  second  vector  and  summing.  In  this  way,  each  vector  is 
weighted  by  its  expected  length,  and  the  vectors  are  progressively  rotated  until  the  largest  vec¬ 
tor  sum  I  is  obtained.  This  occurs  when  each  second  vector  is  the  complex  conjugate  of  the 
first. 

The  weight  function  S(t)  depends  only  upon  the  shape  of  the  spectrum  and  hence  can  be  ob¬ 
tained  from  the  experimental  observations.  In  the  case  of  the  Millstone  Hill  data,  a  mean  of  the 
correlation  functions  c(t)  obtained  in  four  successive  8-minute  runs  is  taken  as  the  weight  func¬ 
tion  employed  in  Eq.  (19).  As  a  trial  value  of  a,  the  value  obtained  from  the  least -mean-squares 
fit  to  the  phase  data  (above)  is  first  employed.  The  sum  I  is  then  recalculated  for  values  of  a 
differing  by  ±5,  ±10,  and  ±15  percent.  Finally,  a  parabola  is  fitted  to  the  seven  values  of  1(a) 
and  the  peak  of  this  parabola  taken  as  the  best  estimate  of  a.  In  the  event  that  the  peak  of  the 
parabola  lies  outside  the  range  ±15  percent  of  the  initial  value,  only  the  least-mcan-squares 
estimate  of  the  drift  is  accepted. 

In  practice,  the  optimum  and  least-mean-squares  straight-line  procedures  are  found  to  yield 
almost  identical  results,  so  that  the  additional  calculations  to  compute  the  optimum  estimate  of 
a  may  be  unwarranted.  Thus,  a  shorter  version  of  the  program  to  extract  the  drift  information 
has  been  written  that  employs  only  the  least -mean-squares  approach.  This  version  has  been  in¬ 
corporated  into  the  program  that  analyzes  the  density  and  temperature  results  as  described  above. 

V.  RESULTS 

A.  Electron  Density 

The  final  electron  density  profile  is  made  available  on  a  high-speed  printout  in  the  form  of 
log.^N  as  a  function  of  altitude.  For  altitudes  between  100  km  and  h  +  100km,  the  A-mode 
results  (which  yield  a  point  every  7.5km)  are  employed  without  any  smoothing.  At  higher  alti¬ 
tudes,  the  smooth  curve  fitted  to  the  combined  A-,  B-,  and  C-mode  data  is  used.  Points  are 
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tabulated  every  10  km  up  to  600km,  and  subsequently  every  25km  up  to  1000-km  altitude.  Any 

results  obtained  for  the  region  above  1000  km  are  tabulated  at  50-km  intervals.  All  these  data 

are  also  written  on  a  "plotting"  tape  and  a  master  data  tape.  The  plotting  tape  is  used  to  drive 

a  CalComp  plotting  machine  to  yield  plots  of  the  logarithm  of  the  density  over  the  height  range 

3  7  3 

100  to  1100  km.  The  density  scale  runs  from  10  to  10  electrons/cm  .  Figures  7(a-f)  provide 
examples  of  plots  produced  in  this  way. 

The  master  data  tape  is  employed  to  accumulate  results  made  over  a  long  period  of  time, 
so  that  it  should  eventually  be  possible  to  search  for  seasonal  variations  and  perform  other  ex¬ 
aminations  of  the  data  that  depend  upon  having  a  statistically  large  sample. 

In  Figs.  7(a-f),  the  presence  of  ground-clutter  echoes  {Sec.  II-B)  at  ranges  near  150  km  has 
given  rise  to  spurious  density  values  between  130  and  160  km.  For  the  daytime  measurements, 
it  is  possible  to  extrapolate  through  this  interval  using  as  a  guide  the  measurments  near  120  km. 
Below  this  altitude,  the  measurements  are  again  suspected  owing  to  the  presence  of  further  clut¬ 
ter  echoes  near  100  km;  these  are  clearly  evident  at  night. 

In  some  instances,  the  second-order  function  fitted  to  the  results  above  h  does  not  de¬ 
max 

crease  monotonically  with  altitude.  This  could  result  from  the  presence  of  spurious  points  caused 
by  satellite  echoes  that  were  not  removed  from  the  data.  In  these  cases,  the  computation  of  the 
density  is  made  to  halt  when  the  gradient  dN/dh  becomes  positive.  In  Figs.  7{d-f),  the  density 
plot  is  not  carried  all  the  way  to  1100km  for  a  different  reason;  namely,  in  order  to  indicate  the 
quality  of  the  results,  the  curve  is  halted  at  the  altitude  at  which  the  uppermost  valid  data  point 
was  obtained. 

The  absolute  accuracy  of  the  density  determination  depends  to  a  large  extent  on  the  reliabil¬ 
ity  with  which  fQF2  has  been  determined.  For  this  reason,  we  plan  to  carry  out  the  analysis  of 
most  of  the  data  employing  the  option  of  substituting  new  values  of  fQF2  {on  punched  cards)  for 
those  recorded  in  real  time.  These  values  can  be  obtained  by  plotting  fQF2  obtained  from  the 
ionosonde  film  together  with  values  obtained  at  nearby  stations,  viz.,  Billerica,  Massachusetts 
{42.6°N,  71.3°W),  Ottawa,  Canada  (45.4°N,  75.9°W),  and  Wallops  Island,  Virginia  (37.9°N, 

7  5.  5°  W).  A  smooth  curve  is  then  drawn  to  follow  the  Millstone  values  except  where  an  obvious 
misinterpretation  has  been  made.  This  procedure  probably  provides  fQF2  accurate  to  0.2  MHz, 
so  that  at  night  the  error  in  the  densities  might  reach  15  percent  from  this  source  alone.  Errors 
in  the  shape  of  the  density  profile  stemming  from  the  scatter  of  the  points  can  be  examined  by 
comparing  successive  profiles.  For  the  most  part,  these  appear  to  be  small,  although  in  some 
cases  a  difference  that  increases  with  altitude  is  noted  between  successive  profiles.  At  1000 -km 
altitude,  these  random  errors  are  thought  to  be  of  comparable  magnitude  to  the  uncertainty  in 
the  absolute  value  of  ^maLX>  Thus,  at  night  it  is  possible  that  the  uncertainty  in  the  density  at 
1000  km  becomes  as  large  as  20  percent. 

B.  Electron  and  Ion  Temperatures 

The  electron  and  ion  temperatures  obtained  for  each  spectrum  are  tabulated  on  the  high¬ 
speed  printout  as  a  function  of  the  delay  of  the  spectrum.  These  same  data  are  written  on  the 
plotting  tape  and  the  master  data  tape. 

When  plotting  the  spectrum  results,  different  symbols  are  employed  to  distinguish  the  0.5- 
and  1.0-msec  results.  A  point  is  also  added  at  120  km  equal  to  the  fixed  temperature  {355°K) 
assumed  in  CIRA  1965  models.  This  serves  to  guide  the  eye  over  the  interval  below  225km.  The 
scale  for  the  temperature  plot  is  chosen  from  a  binary  series  {1000°,  2000°,  4000°,  8000°  K)  in 
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Fig.  7(a-f).  Electron  density  (solid  curve),  electron  temperature  (crosses),  and  ion  temperature  (triangles 
and  diamonds)  obtained  with  new  system  during  24-hour  observing  period  on  9-10  September  1969. 
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Fig.  7(a-f).  Continued. 
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Fig.  7( a-f).  Continued. 
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Fig.  7( a-f).  Continued. 
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Fig.  7(a-f).  Continued. 
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Contin 


order  to  be  able  to  accommodate  the  highest  value  available  for  the  plot.  The  temperature  ratio 
(T0/T.)  results  are  also  plotted  on  a  fixed  scale  running  from  1.0  to  4.0.  Figures  7(a-f)  provide 
examples  of  the  temperature  variation  observed  at  various  parts  of  the  day  for  the  observations 
made  on  9-10  September  1969. 

A  number  of  features  of  these  data  deserve  comment.  Above  ~500  km,  the  electron  temper¬ 
ature  should  chiefly  be  governed  by  conduction  of  heat  from  the  protonosphere.  Thus,  the  tem¬ 
perature  would  be  expected  to  increase  smoothly  with  a  slope  dT^/dh  Up  to  an  altitude 

of  900  km,  most  of  the  measurements  exhibit  this  behavior,  but  Tg  appears  to  increase  in  an 
unexpected  fashion  at  higher  altitudes.  In  order  to  reduce  the  experimental  uncertainty,  we 
have  plotted  in  Figs.  7(a-f)  only  the  average  of  the  temperatures  available  for  97  5,  1050,  and 
112  5km.  Despite  this,  a  systematic  deviation  is  usually  observed.  We  have  not  yet  determined 
the  exact  cause  of  this  phenomenon,  but  believe  that  it  could  stem  from  the  fact  that  Eq.(B)  may 
over-correct  the  temperature  values  as  a  —  1.0  [Eq.  (7)].  In  order  to  circumvent  this  difficulty, 
corrections  to  the  temperature  results  are  not  made  if  ^,0.5,  and  all  such  points  are  then 
discarded.  Thus,  we  would  not  expect  Eq.  (8)  to  introduce  a  large  error.  The  effect  of  an  in¬ 
creasing  abundance  of  light  ions  should  be  to  increase  T.  and  possibly  lower  Te/T^,  but  this  does 
not  seem  to  be  what  is  occurring.  It  suffices  to  say  that  this  is  one  aspect  of  the  reduction  that 
warrants  further  investigation. 

A  second  unsatisfactory  feature  of  the  temperature  results  is  that  at  night  there  is  a  system¬ 
atic  difference  between  the  values  for  T.  derived  in  the  B-  and  C-modes.  This  is  largely  caused 
by  different  estimates  of  T^/T.  for  the  two  pulse  lengths.  In  the  daytime,  on  the  other  hand,  the 
results  agree  remarkably  well.  This  suggests  that  the  difference  is  not  caused  by  instrumental 
effects,  but  is  introduced  in  the  data  analysis.  The  source  of  the  discrepancy  may  be  the  equa¬ 
tions  employed  to  derive  Tp  and  T.  from  the  parameters  (width  and  ratio)  of  the  spectrum. 

These  equations  may  be  inaccurate  if  the  responses  of  the  filters  differ  substantially  from  Eq.  (2). 
The  error  introduced  might  be  expected  to  be  largest  at  night  when  the  width  of  the  spectrum  is 
least.  This  point  deserves  further  attention. 

At  present,  the  accuracy  of  the  temperature  values  seems  limited  by  these  two  systematic 
effects.  During  the  daytime  at  altitudes  in  the  range  225  to  600  km,  the  scatter  in  the  tempera¬ 
ture  measurements  seems  to  be  5  percent  or  less,  suggesting  that  the  random  errors  of  meas¬ 
urement  are  small. 

C.  Vertical  Velocity 

Results  for  the  vertical  velocity  obtained  from  the  Drifts  program  are  tabulated  by  the  high¬ 
speed  line  printer  as  a  function  of  delay,  and  are  written  on  a  plotting  tape.  These  data  are  then 
usually  plotted  as  a  function  of  local  time.  Both  the  least -mean-squares  and  optimum  estimates 
(Sec.IV-D)  can  be  plotted.  Figures  8(a-c)  provide  examples  of  the  drifts  observed  at  three  al¬ 
titudes  during  the  measurements  made  on  25-26  March  1969. 

The  accuracy  of  the  velocity  measurements  has  been  examined  in  detail.*  5  it  appears  that 
the  uncertainty  increases  with  altitude,  being  ~±4 m/sec  at  450km  and  ~±10m/sec  at  900  km 
during  the  day,  and  *6  to  ±2 5 m/sec  over  this  range  at  night.  These  uncertainties  are  set  almost 
entirely  by  the  signal -to-noise  ratio  achieved  in  the  measurements,  and  do  not  depend  greatly 
on  the  variation  with  time  of  the  spectrum  width  or  shape. 
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on  25-26  March  1969.  Crosses  represent  values  obtained  by  least-mean-squares  method,  and  diamonds 
those  from  optimum  method  (Sec.  IV-D). 
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Fig.  8(a-c).  Continued. 
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